Competition between excitonic gap generation and disorder scattering in graphene 
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We study the disorder efTect on the excitonic gap generation caused by strong Coulomb interaction 
in graphene. By solving the self-consistently coupled equations of dynamical fermion gap m and 
disorder scattering rate F, we found a critical line on the plane of interaction strength A and disorder 
strength g. The phase diagram is divided into two regions: in the region with large A and small g, 
771 7^ and r = 0; in the other region, m — and F 7^ for nonzero g. In particular, there is no 
coexistence of finite fermion gap and finite scattering rate. These results imply a strong competition 
A—^ , between excitonic gap generation and disorder scattering. This conclusion does not change when an 

»vj ' additional contact four-fermion interaction is included. For sufficiently large A, the growing disorder 

may drive a quantum phase transition from an excitonic insulator to a metal. 
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I. INTRODUCTION 



The low-energy elementary excitations in undoped graphene are massless Dirac fermions, which are believed to 
1^ ' be responsible for the unusual physical properties of this material ^^. For clean undoped graphene, the density 
$_i of states (DOS) of fermions N{uj) vanishes linearly near the Dirac point, so the Coulomb interaction between Dirac 

c/3 _ fermions is essentially unscreened [ll, [sHl- When the Coulomb interaction is sufficiently strong, A > Ac, it is able to 
(_i ■ open a finite gap for the Dirac fermions by forming stable excitonic pairs. Once this happens, the graphene then has 
2 ' an insulating ground state [y, l8l-[la| . This exotic insulator is induced by the strong electron-electron interaction and 
^ . can be classified as a Mott insulator. 

I ' Although the excitonic insulating transition is a well-studied problem in graphene, so far little work has been 

devoted to the disorder effect on this nonperturbative phenomenon. In any realistic graphene, there are always some 
kinds of disorders, which affect significantly the low-energy transport properties of Dirac fermions. In practice, the 
^ ' scattering of Dirac fermions by disorder can result in important consequences [15| . On the one hand, the disorder 
'— '• scattering induces a finite density of states, N{0), at the Fermi level. The finite N{0) leads to static screening of 
.,j. I the Coulomb interaction. Because the long-range nature of Coulomb interaction plays an essential role in producing 
►^ , excitonic pairing [l5|, the fermion gap decreases with the growing disorder strength g and might finally close when 
QQ ' g is greater than some critical value gc- On the other hand, once the Dirac fermions acquire a finite gap due to 
(T^ , excitonic pairing, the scattering rate and low-energy DOS will be quite different from those of gapless Dirac fermions. 
l/~) ' Therefore, fermion gap generation and disorder scattering are not independent. On the contrary, they have dramatic 
\l ] effects on each other and thus should be treated within a unified formalism. Technically, the fermion gap generation 
[~^ ' can be studied by the Dyson- Schwinger gap equation approach while the disorder scattering rate can be calculated 
^^ ' by means of self-consistent Born approximation (SCBA). In order to examine the disorder effect on excitonic pairing, 
^D I these two approaches should be combined properly. 

'^ ^ ■ In this paper, we examine the disorder effect on the excitonic fermion gap generation caused by Coulomb interaction 

^ , in graphene. By solving the self-consistent equations for dynamical fermion gap m and disorder scattering rate F, we 

found a critical line on the plane of interaction strength A and disorder strength g. This line separates two regions. 

. J In the region with large A and small g, m ^ Q and F = 0, so the ground state is an excitonic insulator with zero dc 

jrt ■ electric conductivity. In the other region, to = and F 7^ for nonzero g, so the Dirac fermions remain gapless and 

" ■ ' acquire a finite scattering rate. In this region, the dc electric conductivity is cr = Ae'^/irh to the leading order of Kubo 

formula, which is exactly quantized and independent of disorder. Apparently, the excitonic fermion gap generation 

due to Coulomb interaction competes strongly with the fermion damping effect caused by disorder scattering. As a 

consequence, the excitonic fermion gap and disorder scattering rate can not be finite simultaneously in graphene. If 

we fix the interaction parameter A at a sufficiently large value, then the graphene may undergo a first-order phase 

transition from an excitonic insulator to a metal at certain critical value gc- These results indicate that the excitonic 

insulating phase transition can take place only in sufficiently clean graphene. 

In addition to Coulomb interaction, there might be contact four-fermion interaction in realistic graphene. Such 
interaction is also able to generate a finite dynamical fermion gap when its strength is sufficiently large. In this paper, 
we also include the four-fermion interaction term in our self-consistent analysis and find that it does not change the 
conclusion. In particular, the excitonic fermion gap generation can not coexist with disorder scattering even when 
additional four-fermion interaction is taken into account. 



The paper is organized as follows. In Sec.2, we study the self-consistent equations of excitonic fermion gap and 
disorder scattering rate when there is only Coulomb interaction. A strong competition between excitonic fermion gap 
generation and disorder scattering is found and a phase diagram on the (A, 5) plane is presented. In Sec. 3, we study 
the self-consistent equations in the presence of an additional contact four-fermion interaction. We show that such 
interaction does not change our phase diagram. In Sec. 4, we discuss the physical implications of our result. 

II. SELF-CONSISTENT ANALYSIS OF EXCITONIC TRANSITION AND DISORDER SCATTERING 

The Hamiltonian of massless Dirac fermion with Coulomb interaction in graphene is given by 

H = VF^ ifja{r)ij -Wipair) 

^ r - A - 

+ irE/ ^.(r)7oV'.(r)T— -^^.'(r')7o^.'(r'), (1) 

a,cr' "''■''"' I I 

where the effective interaction strength is defined as A = 27re^/w_Fe with fermion velocity vf and dielectric constant e. 
As usual, we adopt the four-component spinor field tp to describe Dirac fermion and define the conjugate spinor field 
as -0 = V'Wo- The 4x4 7-matrices satisfy the Clifford algebra [^, Ho ? ]. The fermion flavor is taken to be N and we 
will perform 1/A'^ expansion since Coulomb interaction strength may be too large to be an expansion parameter. The 
disorder potential U{r) couples to fermion as 

N 

Hdis ^Y.'^P C/(r)^.(r)7oV.(r), (2) 

<T=1 "''■ 

where {U{r)) = and (C/(r)[/(r')) — gS{r ~ r'). Here, U{r) plays the role of a random chemical potential and g is the 
disorder strength parameter. For notational convenience, in the following discussion we will assume that vp — 1 and 
restore vp whenever necessary. The total Hamiltonian preserves a continuous chiral symmetry tp —>■ e^^^^ip, where 75 
anticommutes with 7^, which will be dynamically broken once a nonzero fermion gap is generated. 
The free propagator of Dirac fermion is 

Gq ^(iw„, k) = iw„7o - 7 • k, (3) 

where ujn = (2n-|- 1)it/ (i with /? = l/T and n being integers. The Coulomb interaction and disorder scattering modify 
it to the form 

G'^^(za;„,k) = (ia;„ -f isgn(w„)r)7o - 7 • k - m, (4) 

where T is the scattering rate due to interaction. The fermion gap m is generated by the excitonic pairing caused by 
Coulomb interaction and can be obtained from the following Dyson- Schwinger gap equation 

^^7oG(fc)7onP-fc), (5) 

to the leading order of 1/A^ expansion. The unscreened Coulomb interaction is Va{q) — -At. After including the 
dynamical screening effect from collective particle-hole excitations, the effective interaction function has the form 

Vq (9) + n(g) 

The polarization function is defined as 



It has the simple form 



n(g) = -N I j-^Tr [7oG(fc)7oG(fc + q)] . (7) 



N q2 

n{q) = -^^==, (8) 



to the leading order of \/N expansion. It vanishes hnearly as q — >■ in the static hmit q^ — 0, so the Coulomb 
interaction remains long-ranged. For physical flavor N — 2, a. sufficiently strong Coulomb interaction can trigger 
excitonic pairing instability, which generates a dynamical gap for Dirac fermions [0, S fl0l4l3l [la ? ? , [la] ■ 

When disorder scattering is included, a finite fermion damping rate T may be generated, which leads to two 
consequences. First, it changes the energy spectrum of fermions. Second, it yields a finite density of states at Fermi 
energy, which then screens the long-range Coulomb interaction. Both these effects are important and have to be 
incorporated into the gap equation. Because the sign of F is determined by frequency a;„, the gap equation becomes 
much more complicated than that without F. To make our theoretical analysis tractable, here we adopt the frequently 
used instantaneous approximation |8l[l(]|? ]. Taking trace on both sides of Eq.(5), we get 

1 ■r-^ /■ rf^k m(iu;„,k) 



m(ie„,p) = ^J2 



(27r)2 (iiOn + irsgn(a;„))2 + jkp + m2(iw„, k) 
xF(ie„-iw„,p-k,F). (9) 

It is possible to sum over w„ only when the frequency-dependence of fermion mass function is neglected, which is the 
so-called instantaneous approximation '8|, llOl llSl ? ] . In this approximation, this equation becomes 



-(p) - -pUJ 



d^k m(k) 



(27r)2 {iuj„ + iFsgnK))2 + |k|2 + m2(k) 
xF(p-k,F). (10) 



After carrying out the frequency summation, we have 

d^k m(k) 



m(p) 



(27r)2 ^|k|2+m2(k)7r 



2 27rT 27rT 

To make an analytical computation, we take the limit T — > 0, so that 



I F V|kp+m2(k)\ , , 

xI"i^U + 7r:77; + ^ ^' '._^ ^M y(p-k,F). (11) 



,. T ,/l r . v/k2 + m2 k \ / J\i^ + m? k \ 

lim Imib - H h i- —^ = arctan -^ ^-^ , 12) 

T^o ^ I 2 27rT 2ttT \ F i ' ^ ' 



which then leads to 



, d^k m(k) 1 /V|k|2-^m2(k) 

fnin) ~ / — — arctan ' ^ 

' (27r)2^|k|2+„i2(k)^ 



1 



^ + n(p-k,F)' 



(13) 



We now need an analytical expression for the polarization function n(q, F), which is crucial to describe the screening of 
long-range Coulomb interaction from particle-hole excitations. In the Matsubara formalism, the polarization function 
within instantaneous approximation is given by 



+00 ,2 



nj TJ^ ^ ,2j^ 
n(q,r,F)=-- J2 / 7^^Tr[Go(z^„,k)7oGo(j^„,k + q)7o]. (14) 



n— — oo 



After introducing Feynman parameter and summing up imaginary frequencies, the polarization function is found to 
have the form 



2N I I C 

n(q,r,F) = ^- / dx dt<-flm 

TT Jo JC„ * 



'1 F-fit^ 
^^2 + 1^^ 

i^-G2aim[^(i + iif)]l 



8t ( ' ^^^) 



where Cq — \/x{l — a;)|q| and an ultraviolet cutoff A is introduced. This expression is very complicated, but can be 
further simplified in the zero temperature limit. As T — ?> 0, we have limr^o Im [V-" (5 + 3^)] — arctan (1?). After 
tedious but straightforward computation, we finally have 



n(q,r) = 



2N 



rin 



F(|q|,r) + |q| / dx^JxiX - x) 
Jo 



X arctan 



'y/x{l~x)\q\ 



(16) 



where the function F{x) is 



^(|q|,r) 



-r 



rV4r^-|qP 



arctan 



y4r^-|q| 

_^ _ rV|qP-4n j^ / |q|-v/|qP-4r2 
2|q| \^|q|+^|q|2-4r2 

-r. 



o< ^ <2 

^>2 

|q| _o 
r - ^ 



(17) 



If we further let |q| = 0, the polarization function takes a finite value 

n(o,r) 



^rm^^ 



(18) 

which implies that the long-range Coulomb interaction now becomes short-ranged. In fact, 11(0, F) is precisely the 
Dirac fermion DOS at zero energy iV(0) induced by disorder scattering [l3|. As shown previously [l^, the long-range 
nature of Coulomb interaction plays an important role in driving excitonic gap generation. Once the long-range 
Coulomb interaction is statically screened by disorder scattering, excitonic gap generation may be suppressed. 

In the above equations, the scattering rate T is simply assumed to be an arbitrary constant. However, T reflects the 
fermion damping due to disorder scattering and may depends heavily on the fermion gap. Indeed, fermion scattering 
rate F and fermion gap m{p) can affect each other. In a more refined treatment, F should be calculated on the same 
footing with fermion gap m(p). 

One frequently used approach of calculating F is to first average over random potential U(r) and then make mean- 
field analysis in the saddle point approximation, which amounts to SCBA J18l - l2]| . Although SCBA was questioned 
recently [22| , at present there is no better choice when dealing with the interplay of disorder and excitonic pairing. 
Within this approximation, the disorder scattering rate is determined by the following equation 



F(w„) 



-Tr 



70 



d^k 



(27r)2 {iujn + isgn(a;„)F(w„))7o - 7 • k - m{k) 



(19) 



In principle, the scattering rate F receives contributions from both elastic scattering by quenched disorder and inelastic 
scattering by Coulomb interaction. However, the later contribution is forced to vanish at zero energy by the phase 
space restriction originated from the Pauli exclusion principle. However, the scattering rate from disorder can be 
finite even at zero energy. Thus, we simply ignore the Coulomb interaction in this equation. 

From the above equation, we know that F(a;„) does not depend on momentum k because any k-dependence will be 
lost after integration over k. Further, no energy is transferred during the disorder scattering process, so the scattering 
rates F(w„) at different w„ are completely independent. Thus, we can drop the energy dependence of F(a;„) and focus 
only on the equation for F(0). In fact, in the derivation of gap equation, we have already assumed that the scattering 
rate F is independent of both momentum and energy (otherwise it is impossible to get an analytical expression for gap 
equation). Now we see that this assumption is reasonable. Taking advantages of these simplifications, we eventually 
have 



9 



'-tJ * 



T2 



m^Ck) 



(20) 



When m — 0, the integration over k can be precisely done yielding a scattering rate F — Ae^^'^/^, which has been 
known for many years [l^. If we require that the function m(k) be given by fermion gap equation Eq.(13), then we 
obtain two self-consistently coupled equations for dynamical fermion gap and disorder scattering rate. 

Before doing numerical computation, it is helpful to first make some qualitative analysis. For simplicity, we assume 
a constant fermion itiq. Now the integration over k in Eq.(20) can be performed, yielding the solution 



p2 



A^ 

o4ir/g _ I 



(21) 



From this result, we know that T can have a nonzero solution only when g is greater than some critical value 

47r 
In(rT^) 



9 > , ., . A2. - (22) 



Apparently, only sufficiently strong disorder can induce a finite scattering rate when fermions are massive. From this 
qualitative analysis, we expect a competition between excitonic gap generation and disorder scattering. 



FIG. 1: Phase diagram of graphene on the plane of (A, (;). In region 1, m ^ and F = 0. In region II, jn = and F 7^ when 
(? 7^ 0. This phase diagram is valid for weak disorder. 

In reality, the fermion gap 7TT,(k) is not a constant, but is a function of fermion momentum. It is therefore essential 
to solve the equations (13) and (20) self-consistently. After carrying out extensive numerical computation at physical 
flavor A^ = 2, we found a critical line on the plane of {X,g) and show the phase diagram in Fig.l. In region I with 
large A and small g, the dynamical fermion gap has a nontrivial solution, m ^ 0, but the disorder scattering rate 
have only trivial solution, F = 0. In region II, the fermion gap has only trivial solution, to = 0, while the disorder 
scattering rate acquires a finite value, T ^ (certainly, F vanishes when g = 0). We did not find any evidence for the 
coexistence of nontrivial m and nontrivial F in numerical computation. The numerical results confirm our expectation 
that there is a competition between the generation of excitonic fermion gap and the generation of disorder scattering 
rate: fermion damping F suppresses the formation of excitonic pairs, and fermion gap to, prevents the appearance of 
fermion damping. 

In region I on Fig.l, the excitonic gap generation wins the competition, hence the ground state is an excitonic 
insulator with zero dc electric conductivity u = 0. In region II, the disorder scattering wins, then the ground state 
is gapless and the Dirac fermions acquire a finite scattering rate. In region II, the accurate calculation of electric 
conductivity is a problem in debate [ll, 01 • To the leading order of the Kubo formula, the dc electric conductivity is 
known to be cr = Ae'^/nh, which is independent of disorder and universal [17h21| |. When disorder strength g grows 
from I to II, a jumps suddenly from zero to a universal value at certain critical value gc- This is a first-order excitonic 
insulator to metal phase transition driven by growing disorder. Note, however, that the excitonic phase transition is 
of first-order only in the presence of finite disorder, g ^ 0. In the clean limit, g = 0, the excitonic transition is not 
of first-order but is continuous as the Coulomb interaction parameter A is varied. There is indeed a debate on the 
nature of this transition in the clean graphene. An infinite-order transition was claimed in some theoretical analysis 
[lot Il6| , whereas a conventional second-order transition was found in recent numerical work [l3| . 

As mentioned earlier, the disorder considered in this work is random chemical potential, which may be generated 
by local defects, neutral impurity atoms, or neutral adsorbed atoms in the plane of graphene [j,|j|- When this type of 
disorder is smooth at the atomic scale, it will not mix the two inequivalent Dirac points. The parameter g is indeed the 
product of the concentration of impurity atom (or defect) and the strength of a single impurity atom. In practice, the 
magnitude of the critical disorder parameter gc depends on the fermion fiavor and the Coulomb interaction strength A. 
If we take physical fiavor N — 2 and assume that graphene is suspended in vacuum so that e'^/vp^ ~ 2.16 [23] . namely 
A ~ 13.57, then it is easy to obtain gc ~ 0.36. When graphene is placed on certain substrate, the Coulomb interaction 
strength is reduced by the screening due to substrate. From Fig. 1, we know that critical value gc decreases as A is 
lowered. 

We would like to emphasize the importance of making a self-consistent analysis in our problem. In fact, if we solve 
the fermion gap function (13) by assuming a finite constant scattering rate F, then there is a critical value Tc- The 
gap equation has no nontrivial solution for F > Tc, but a finite fermion gap is opened for F < Fc. In this case, there 
is coexistence of finite fermion gap and finite scattering rate when F < Fc. Similarly, if we solve the SCBA equation 
(20) by assuming a free constant gap rriQ, there will be coexistence of finite fermion gap and finite scattering rate 



when the inequahty (22) is satisfied. In both these cases, there will be the third region with finite fermion gap m 
and disorder scattering rate T lying between region I and region II on the phase diagram. In this third region, the dc 
electric conductivity would be 

TTAl I ^ + TO^ 

which depends on disorder strength and displays metallic behavior even when the fermions are gapped. However, 
when the equations (13) and (20) are solved self-consistently, the fermion gap m and scattering rate F can not be 
finite simultaneously, and can not be zero simultaneously when g ^ 0. As a consequence, the dc electric conductivity 
is either zero or exactly quantized, and does not explicitly depend on disorder strength. 

III. EFFECT OF ADDITIONAL FOUR-FERMION INTERACTION 



In additional to Coulomb interaction, the contact four-fermion interaction may also be important in realistic 
graphene. In the language of field theory, such contact interaction is normally described by the Gross-Neveu model 
[2J, [23. It can make additional contributions to the generation of finite excitonic fermion gap 24, 25|]. A natural 
question is whether the contact interaction alters the phase diagram shown in Fig.l. When a four-fermion interaction 
is included in our self-consistent analysis, it is in principle possible to obtain a coexistence of finite excitonic fermion 
gap and finite scattering rate. In this section, we will examine this possibility. 

As an concrete example, we consider the following Gross-Neveu model 

r ^ r - 
Hgn ^ J-J2 {Mr)Mr))' ■ (24) 

This interaction term does not respect the continuous chiral symmetry, but respects a discrete chiral symmetry. 
Its role in excitonic gap generation was analyzed in recent years |lll . Il2l llSl Il6l | . Unlike Coulomb interaction, the 
interaction strength does not depend on fermion momentum and energy and there is no dynamical or static screening. 
We first ignore the Coulomb interaction and consider the Gross-Neveu model only. The corresponding gap equation 
is 



77l(p) 




(25) 



where g' = GNA. From this equation, it is clear that fermion gap is indeed independent of momentum. At zero 
temperature, we can reduce gap equation to 



1 . ,^ r«.i^,^^.„ta„^^Et±i! , „o, 



In the clean limit, F = 0, then 



1 = T^T-^ / rf|k||k| , ^ ^, (27) 

AAT^X I I II I /n |2 . 9' ^ '' 



which has solution 



f 4NTT ^ 



m 

g' 



(28) 



There is a critical coupling g'^ — ANtt. When g' < g'^, there is no non-trial solution for m; when g' > g'^, there is 
non-trivial solution for m. In the presence of disorder, the scattering rate F is given by the SCBA equation (20). 
Since now fermion gap tti is a constant, it is easy to obtain the following expression 

r' = ^r4- m^- (29) 



After solving fermion gap equation (26) using this T, we found no coexistence of finite fermion gap and finite scattering 
rate. 

We next consider the case when both Coulomb and Gross-Neveu interactions are present in graphene. Recently, 
Gamayun et al. showed that the analytical results can agree with numerical simulation results once the Gross-Neveu 
model is included [16i] . Now the whole gap equation has the form 



m(p) = / 



d^k m(k) 1 

47r2 ^|k|2 + n^ ip_M + n(o, p - k, r 



xiarctan vWTZM 

TT \ r 



^' f'^ ™('') =iarctan^S!+ZM , (30) 



A^A J 47r2 ^|k|2+m2(k) tt ^ T 

which couples self-consistently to SCBA equation (20). After numerically solving these equations, we found no 
coexistence of finite fermion gap and finite scattering rate. Although the magnitude of fermion gap in the excitonic 
insulating phase becomes larger after Gross-Neveu interaction is included, the qualitative phase diagram shown in 
Fig.l does not change. 

IV. SUMMARY AND DISCUSSION 

In this paper, we studied the disorder effect on excitonic fermion gap generation, i.e., excitonic insulating transi- 
tion, due to Coulomb interaction in graphene. By solving the self-consistent equations of fermion gap and disorder 
scattering rate, we found a strong competition between excitonic fermion gap generation and disorder scattering. As 
a consequence of this competition, fermion gap generation can not occur simultaneously with disorder scattering in 
graphene. The phase diagram presented in Fig.(l) is the main new output of this paper. We also showed that this 
phase diagram is not changed by additional contact four-fermion interaction. 

In any realistic graphene, the Dirac fermions are always scattered by some kinds of disorder. Our results indicate 
that, even if the Coulomb interaction is indeed sufficiently strong, the excitonic transition will be completely suppressed 
if the disorder strength is not small enough. Apparently, the excitonic transition can most possibly be observed in 
very clean graphene. 

We would point out that our results are obtained using SCBA. It would be interesting to examine the effect 
of higher order corrections. The most important corrections come from the fluctuation effect associated with the 
Anderson localization. When F ^ 0, the metallic electric conductivity a — Ae^/nh is subjected to diffusion and 
Cooperon vertex corrections [261 ] , which will lead to Anderson metal- insulator transition. When Coulomb interaction 
is absent, A = 0, the localization of massless Dirac fermions was found |2j,[23|. For small but finite A, though Coulomb 
interaction is not strong enough to open fermion gap, it can have important infiuence on transport properties [y, |28| - 
1341 . In particular, it may destroy the fermion localization [35|. However, there is currently no widely accepted theory 
for the interaction effect on localization [35|. On the experimental side, earlier measurements suggested that the 
undoped graphene exhibits a universal minimum conductivity [i[, ct = Ae^/h. More recently, it becomes clear from 
extensive transport experiments that the minimum conductivity in undoped graphene is not universal but instead 
strongly sample-dependent [4j, i36n38i] . Regardless of the precise value of minimum conductivity, it seems that the 
undoped, gapless graphene is metallic and free of localization at experimentally accessible temperature 0-01 : at least 
for weak disorder. 

In this paper, we are mainly interested in the regime of strong Coulomb interaction. For large A, the localization 
effect is even more involved because it is entangled with the non-perturbative phenomenon of excitonic gap generation. 
As already explained, a self-consistent treatment is crucial in our problem, hence excitonic gap generation and fermion 
localization can not be studied separately. Technically, the fermion gap function m{p) appearing in (13) should be 
used when computing the diffusion and Cooperon vertex corrections, while these vertex corrections should be included 
in the polarization function n(p) and the fermion gap equation (13), which correspond to the Altshuler-Aronov type 
corrections [20| . Unfortunately, although in principle it is possible to examine the importance of Anderson localization 
effect, the self-consistent equations obtained in this way are too complicated to be analyzed theoretically or numerically. 

In this paper, we considered only one particular type of disorder, i.e., random chemical potential. Our self-consistent 
analysis may be extended to study the effects of other types of disorders [331 , such as random gauge field or random 
mass, on excitonic pairing formation. Specifically, ripples are believed by many people to be important in graphene 
and thus attracted intensive investigation ,4^| . Such ripple configuration can be described by a random gauge potential 



8 

|40| . It is interesting to study the ripple effect and to examine whether ripples can drive an analogous phase transition 
in the future. 

We would like to thank P. Pyatkovskiy for pointing out a mistake in our previous calculation of polarization function 
and W. Li and J. Wang for helpful discussion. G.Z.L. is supported by the National Science Foundation of China under 
Grant No. 11074234 and the project sponsored by the Overseas Academic Training Funds of University of Science and 
Technology of China. 
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